arXiv:1508.01257vl [astro-ph.GA] 6 Aug 2015 


Variable Reddening and Broad Absorption Lines in the 
Narrow-line Seyfert 1 Galaxy WPVS 007: an Origin in the Torus^ 

Karen M. Leighly 

Homer L. Dodge Department of Physics and Astronomy, The University of Oklahoma, 440 

W. Brooks St., Norman, OK 73019 

Erin Cooper 

Homer L. Dodge Department of Physics and Astronomy, The University of Oklahoma, 440 

W. Brooks St., Norman, OK 73019 

Dirk Grupe 

Department of Earth and Space Science, Morehead State University, 235 Martindale Drive, 

Morehead, KY 40351 

Donald M. Terndrup 

Department of Astronomy, The Ohio State University, 140 W. 18th Ave., Columbus, OH 

43210 

S. Komossa 

Max-Planck Institut fiir Radioastronmie, Auf dem Hiigel 69, 53121, Bonn, Germany 
Received_; accepted_ 


Accepted for publication in ApJ Letters 




- 2 - 

ABSTRACT 

We report the discovery of an occultation event in the low-luminosity narrow- 
line Seyfert 1 galaxy WPVS 007 in 2015 February and March. In concert with 
longer timescale variability, these observations place strong constraints on the 
nature and location of the absorbing material. Swift monitoring has revealed 
a secular decrease since ~ 2010 accompanied by flattening of the optical and 
UV photometry that suggests variable reddening. Analysis of four Hubble Space 
Telescope COS observations since 2010, including a Director’s Discretionary time 
observation during the occultation, shows that the broad-absorption-line velocity 
offset and the C IV emission-line width both decrease as the reddening increases. 
The occultation dynamical timescale, the BAL variability dynamical timescale, 
and the density of the BAL gas show that both the reddening material and 
the broad-absorption-line gas are consistent with an origin in the torus. These 
observations can be explained by a scenario in which the torus is clumpy with 
variable scale height, and the BAL gas is blown from the torus material like 
spray from the crest of a wave. As the obscuring material passes into our line of 
sight, we alternately see high-velocity broad absorption lines and a clear view to 
the central engine, or low-velocity broad absorption lines and strong reddening. 
WPVS 007 has a small black hole mass, and correspondingly short timescales, 
and so we may be observing behavior that is common in BALQSOs, but is not 
typically observable. 

Subject headings: quasars: absorption lines — quasars: individual (WPVS 007) 
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1. Introduction 


Many quasars exhibit broad (width An > 2000 kms ^), blue-shifted absorption 
lines in their rest-fr ame UV spectra that are evidence for winds emerging from the 


central engine (e.g., 


Wevmann et ah 


19911); these objects are called broad absorption-line 


quasars (BALQs). These winds are an essential part of quasars: they carry away angular 


momentum and thus facili 

the intergalactic medium ( 

tate accretion. 

:hev distribute chemically enriched gas through 

Gavaliere et al. 

2002 

), and they may iniect 

cinetic energy into 

)• 

the host galaxy, influencing its evolution (e.g.. 

Scannanieco & Oh 

2004 


Despite the ubiquity of quasar winds, our understanding of the mass outflow rate, 
geometry, and driving mechanism is limited. Line variability provides a powerful tool 


for studying the absorbing gas. Modest variabili ty on time scales o 


common, seen in ab out one-third of BALQs (e.g.. 


Gibson et al. 


2008 


a year or less is 


Filiz Ak et al. 


Capellupo et al. 


2011 


20131) . Variability is typically limited to changes in line equivalent width. 


sometimes in a single component and often at higher velo cities, but occasiona l 


appearances or disappear ances of components occur (e.g.. 


2011 


Filiz Ak et al. 


20121 ) 


Leighlv et al. 


2009 


V dramati c 


Hall et al. 


Another type of absorption phenomenon, occultation of the X-ray emitting region, 
has recently been observed. One example is NGC 1365, which showed occultation in the 
broad-line region and transitioned between reflection-dominated and Gompton-thin X-ray 


spectra in four days, indicating a compact X-ray emission region flRisaliti et al. 


20071) : many 


^Based on observations made with the NASA/ESA Hubble Space Telescope, which is 
operated by the Association of Universities for Research in Astronomy, Inc., under NASA 
contract NAS 5-26555. These observations are associated with programs 11733, 13015, and 


14058. 
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other examples are known (for a review, iBianchi et al.ll2012l) . Thus, it appears that active 
galactic nucleus (AGN) obscuration is more complicated than the simple unihed model for 
the Seyfert 1/2 dichotomy. 


WPVS 007 i s a low luminosity 


Seyfert 1 galaxy (iGrupe et ah 


[Mv ~ —19.7), nearby {z = 0.02882) narrow-line 


19951) with peculiar UV properties. In a 1996 Hubble Space 


Telescope FOS spectrum, it showed moderately narrow absorption lines that classihed 
it as a mini-BAL quasar. A FUSE spectrum from 2003 showed evidence for a dramatic 
change in the absorption: in additi on to the mini-BAL s, a broad absorption line with 
Gmax ~ 6,000kms“^ had appeared flLeighlv et al.ll2009l) . 


We present Swift UV photometric monitoring, which reveals the hrst UV occultation 
event ever observed. We also discuss four HST GOS spectroscopic observations, including 
one during the occultation, that indicate a correlation between the BAL outflow properties 
and the reddening. Analysis of timescales suggests that the variability has a common origin 
in the torus, and we suggest a scenario to explain the observations. Full details of the 
observations, reduction, and analysis will follow in E. Gooper et al. (2015 in preparation). 


2. Observations 


2.1. Swift Observations 


WPVS 007 has been monitored w ith Swift since 2005 


with the UV-optical telescop e (UVOT, 


near UV and optical bands. 


Roming et al. 


Grupe et al. 


(iGrnpe et al. 


2007 


2008 


20131) 


2005h using its six hlters spanning the 


(l2013f) reported observations through 2013 July 


13. Monitoring continued using UVM2 (rest effective wavelength Aeg = 2159A) through 
the remainder of 2013 and 2014. A secular decrease since ~ 2010 is seen Fig. [H and the 
UV continuum flux dropped signihcantly following the observation on 2014 November 29. 
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A six-filter observation on 2015 Febrnary 28 showed that the spectrnm had become mnch 
redder, so we obtained observations with approximately weekly cadence to the present time 
using all hlters. 


We used the six-point photometry to quantify the spect ral variab ility. The 147 spectra 


were £t simultaneously to a power law plus SMC reddening flPei 


1992I) . with the power law 


index held constant. The results are shown in Fig. [T] We estimate that the continuum was 
occulted for about 60 days based on the magnitude and E[B — V) variability (Fig. [1] inset). 
Color-magnitude plots shown in Fig. [U display the spectral variability, and superimposed 
lines show that the color and brightness changes are consistent with variable reddening. 

Intrinsic optical/UV spectral variability is not a favored explanation for the observed 
behavior. Principally, the recent event has a characteristic occultation p rofile, as seen i n 


other astronomical objects such as Cataclysmic Variables. Additionally, 


Kokubo et al. 


(120141) show that the typical intrinsic-variability difference spectrum is a power law, but 


that model does not £t the color/magnitude variability as well as SMC reddening (Fig. d]). 


2.2. HST Observations 

WPVS 007 has been observed using HST hve times; the four more recent observations 
are marked on Fig. [H Spectra in 2010 and 2013 sample the secular decrease, and the most 
recent observation in 2015 occurred during the occultation. We display the HST spectra 
and averages of the Swift photometry near the HST observations in Fig. [2] (left). 

Fig. [2] (right) shows the C IV spectral region, consisting of the deep mini-BAL and the 
higher-velocity broad absorption lines. The mini-BALs, also present in Si IV, Lya, N V, 
and S III, do not vary, suggesting that they originate in a separate large-scale component, 
unassociated with the BALs; a discussion of their properties is deferred to Cooper et al. 
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(2015, in preparation). 

The spectra reveal several model-independent properties. The broad absorption lines 
show dramatic variability. The C IV BAL depth is never lower than the continnnm level of 
the 1996 spectrum, and the bottoms of the 2013 spectra are flat. This snggests the presence 
of a non-variable continnnm component that is not absorbed by the BALs. The origin of 
this is not clear; some of it is likely to be intrinsic line emission. 


The BAL velocity offsets appear to be correlated with the flux. The variability is 
most clearly observed in the 2010 high-state spectrum, which shows higher 14iax and 14iin 
compared with the others. There is also evidence that the low-state spectrum has weaker 
high-velocity absorption: the 2013-minus-2015 difference spectrum (Fig. [3]) shows the Lya 
line in 2015 has a red wing not present in the 2013 spectra. This suggests that the Lya line 
was absorbed by the high-velocity component of the N V broad absorption line in 2013, but 
that component is weak or absent in the 2015 spectrum. 


The C IV emission-line width decreases as the reddening increases (Fig. [3]). The 
difference between earlier spectra and the low-state spectrum shows positive residuals 
bracketing the rest wavelength. This is most clearly seen in the 2010 spectrum, where 
the large BAL leaves the blue wing of the emission line intact, but it can also be 
seen in the red wing in the 2013 spectra. Spectral htting conhrms the FWHM variability: 
FWHM(C IV)= 3275 kms“^ in 2010, ~ 1850 kms“^ in 2013, and = 1450 kms“^ in 
2015. This suggests that the same reddening/occultation causing the continuum spectral 
variability also covers the inner broad-line region. 


All of the HST spectra were ht simultaneously using Sherpa (iFreeman et ah 


200lh . 


There is structure in the continuum that is most prominent in the 2015 spectrum where the 
AGN continuum is most absorbed. Since the 2013 minus 2015 difference spectra are smooth, 
we suggest that this structure is weak line emission; there is moderate Fe II emission 
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at longer wavelengths. We created a constant continuum component by modeling the 
2015 spectrum empirically with broad features, subtracting the model from the spectrum, 
then htting the residuals with a high-order spline, and adding the model to a constant; 
see Cooper et ah in prep, for details. The AGN central engine continuum was modeled 
with a broken power law. The emission lines were modeled using Lorentzian prohles. 

The mini-BALs were modeled using four sets of gaussian optical-depth absorption lines, 
and assumed to absorb emission lines and both continua. The BALs were modeled using 
gaussian optical-depth absorption lines absorbing the AGN continuum and emission lines, 
but not the constant continuum component. The BAL offsets and widths were tied together 
among the C IV, Si IV, and N V lines for the 2013 and 2015 spectra. The resulting £t is 
shown in Fig. [3l 


3. Discussion 

Our new observations of WPVS 007 confirm its unique nature among AGN. Swift 
observations of short-timescale optical/UV spectral changes reveal the hrst occultation 
event observed in the UV. Long timescale reddening changes appear to be correlated with 
unusual changes in the broad absorption-line velocities and the broad emission-line widths 
observed by HST . In this section, we argue that the occulting material and absorbing 
gas have their origin in the torus, and present a scenario that can explain the variability 
behavior. 


3.1. Size Scales in WPVS 007 


Fig. m s 
hole mass. 


rows radii an d size scales for WPVS 007 as a function of plausible black 
(2009) estimated the black hole mass is 4.1 x 10® Mq. 


Leighlv et ah 










Vestergaard fc PetersonI (120061) estimated that masses derived from scaling relations 


are nncertain to 0. 


Leighlv et al. 


3 dex, so we consider masses between 1.5 x 10® and 1.1 x 10^ Mq. 


(120091) performed photoionization analysis of the FUSE spectrnm, using an 
X-ray-weak spectral energy distribution with anx = ~l-9, i nferred from the hard X-ray 


detection in a long observation by Swift (iGrnpe et al 


20081 ). Scaling this SED with the 


1996 HST spectrum indicates a bolometric luminosity of 5.2 x 10^® ergs 


We use a sum-of-blackbodies accretion disk with rj = 0.1 (e.g., iFrank et ahll2002l) to 
estimate the radii characteristic of continuum emission at 1550A (i?i 55 o, under C IV), 
2150A (i? 2 i 505 rest-frame Agg for the UVM2 filter), and 3350A (i? 335 o, rest-frame Agfj for the 
U filter). 


Leighlv et al.l (1200911 estimated that the H/3 emission 


region size 

Rna) 

is 

9.7 X 10 ® pc 

Bentz et al. 

(2006 

)• 

Grupe et al. 


based on the 1996 HST spectrum and the relations from 

ri 

(1201311 foun d a sm aller value of 4.2 x 10 ^ pc based on the Swift photometry and using 


Kaspi et al. 


(| 2000 |). 


The H(5 line width is 1190 km s ^ (iLeighlv et al. 


2009 1. The C IV line width varies. 


becoming narrower in the low state, which we attribute to reddening of the inner broad-line 
region by the same material responsible for the long-term secular decrease. To estimate 
the inner BLR radius RciVi we simultaneously ht the low-state spectrum with a single 
Lorentzian prohle, and the other three with broad and narrow Lorentzian components. The 
width of the broad component is found to be 4170 kms“^. We estimate Rciv to be between 
3.4 X 10“^ and 7.9 x 10“^ pc by scaling with Rhp and assuming that the BLR is virialized. 

The torus is estimated to he between the dust sublimation radius and the radius of its 


peak emission at 12yum. T 


le dust sublimation radius can be estimated to be the 


reverberation radius Rtk ( Kishimoto et al. 


<-band 


2011 1. Using Eq. 1 of 


Kishimoto et al 


(1201 ih 


and the 5500A flux density from the 1996 HST spectrum yields Rtk = 0.036 pc. The 
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12fim half-light radius is best measured by interferom etry; however, /‘),(1 2u'm) does not 


have a clear luminosity scaling relationship like Rtk flBurtscher et ah 


2 OI 3 I). We estimate 


a plausi ble range for 


those in 


Burtscher et ah 


12jjm) by comparing the WPVS 007 bolometric luminosity with 


fl2013[) Table 6, estimating that Ri/2{12fim) may be between 0.3 


and 2 pc. The ratio of our estimated R^ / 9. (12um) and R tk ranges from 8 to 56, consistent 


with the value of ~ 30 generally observed flNetzer 


2015 


, their Fig. 4). 


The BAL absorbing gas radius is estimated using the S IVA1063 and S IV*A1073 


broad absorption lines from the FUSE spectrum flLeighlv et al 


2009 ). Knowing that 


the excited-state population increases with density, we follow the standard method (e.g.. 


the excitec 

-state 

Arav et al. 

2013 


Borguet et al. 


20131) and hnd the gas density lies between log(n) = 4.9 


and 5.1 [cm depending on whether partial covering is accounted for (see Cooper et al. 
2015 in preparation for details). The absorber distance de pends on the photoionization 


parameter through U = Q/dii R'^n^.c; 


Leighlv et al. 


(120091) showed log 7/ —0.3. The 


upper limit cannot be determined with the spectra in hand; very hig. 
chara cteristic of log 17 > 0 he in the extreme UV (e.g.. 


liand; very hig 

1 ionization lines 

Hamann et al. 

199§; 

Arav et al. 


20131). However, the presence of P V in the FUSE spectrum indicates a high column density 


(e.g., 


Hamann 


19981) . A higher ionization parameter would drive the zone deeper into 


the gas, increasing the inferred Nh and associated kinetic luminosity. Assuming that the 
ratio of the kinetic luminosity to bolo metric luminosity is implausibly large if > 20%, and a 


global covering fraction 72 = 0.2 (e.g.. 


Dunn et al. 


20101) ■ we estimate a rough upper limit 


on the ionization parameter. These constraints yield Real between 0.17 and 1.47 parsec, 
consistent with the torus. 
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3.2. The Occulting Cloud Location is in the Torus 

We estimate the distance of the occulting cloud from the central engine using the 
occultation timescale of T = 60 days f ^2.1l) . V = R\/T, and the Keplerian velocity. We 
approximate the occulting-cloud size to be the radius of the continuum emission region 
(i.e., the cloud does not fully cover the source), since we observe residual continuum in the 
low-state spectrum. We perform the calculation for the three continuum-emitting radii 
computed in §O(i?i 550 , ^ 2150 , ^ 3350 ), yielding RAbs{i 550 ), i?Abs( 2 i 50 ), and RAbs{ 3350 ) (Fig.S]). 
The occultation is characterized by changes in reddening, and for dust to be present, the 
occulting cloud should lie at i? > Rtk- This criterion is obeyed for i?Abs( 2 i 5 o) and RAbs{i 55 o), 
but not RAbs{ 3350 ), and therefore the occulting cloud can largely cover the inner continuum 
emission region, but not the outer. It seems plausible that the occulting material is a clump 
embedded in a larger region of dusty gas responsible for the longer timescale reddening 
changes, and both are associated with the torus. 


3.3. The Broad Absorption Line Gas Originates in the Torus 

To fully understand BAL variability, we ideally need frequent HST spectroscopic 
monitoring. While we do not have that, our observations, although sparse, were fortuitously 
placed, sampling the long-term secular decrease and increase in reddening between 2010 
and 2015. The apparent correlation between BAL velocity offsets, reddening, and C IV 
emission-line width suggest they all have a common origin. While we have already estimated 
Real, based on the absorbing gas density, and it is consistent with an origin in the torus, 
we also estimate a dynamical radius for inner-BLR occultation Rabs(civ) for consistency. 
A very rough estimate is obtained if we assume that the BAL variability timescale is the 
same as the C IV line-width variability timescale, and require that the BAL absorber 
cover Rciv- The timescale should be longer than the separation between the two 2013 
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observations (190 days), because we see little difference between those spectra, but no 
longer than the separation between the 2013 December observation and the 2015 March 
observation (473 days). The derived range of radii Rabs{civ) using the timescale of 473 
days is also consistent with an origin in the torus, as well as Real (Fig. Hj). 


3.4. A Scenario for the Spectral Variability in WPVS 007 


Typical BAL variability has been attributed to several causes (e.g., ICapellupo et al. 


20121) . Coordinated var iability in multiple t roughs, for example, can be attributed to 


ionization changes (e.g.. 


Hamannetah 


20111) . We have seen, however, that the variability 


in emission-line width, BAL velocity, and continuum reddening in WPVS 007 are unlike 
that seen in other AGN. Moreover, we have shown that the absorption and the reddening 
are plausibly associated with the torus. It seems likely that all of these unusual features 
have a common cause, and we suggest a scenario that can explain them. 


We suppose t hat the torus is not uniform but clumpy, a picture now widely accepted 


(e.g., 


Netzer 


2015 


, for a review), and that the torus scale height is not constant with 


azimuth (Fig. [5]). If our view to the central engine in WPVS 007 skims the edge of the 
torus, rotation of the (clumpy) variable-scale-height material can naturally explain the 
short and long timescale reddening variability. 


We suppose further th at the torus may be the source of the BAL gas. This is not a new 


idea; 


Krolik &: KrissI (120011) suggest an origin of warm absorbers from evaporation of the 


inner edge of the torus. Perhaps the dusty gas is ablated from the torus by the continuum 
radiation like wind ablates spray from the crest of a wave. Dust has a high absorption 
opacity that may help accelerate the gas, and we might expect the outflow velocity to be 
larger farther above the torus (Fig. [5]). When a low-scale-height region orbits into our line 
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of sight, we see through higher-velocity gas and less reddening. When a high-scale-height 
place orbits into our line of sight, we see more reddening and a lower-velocity outflow. This 
scenario also explains the emission-line-width variability if the absorber also covers the 
inner broad-line region. 


4. Summary 


Recent improvement in monitoring capabilities reveals that occultation events are 


not uncommon in the X-ray band (e.g.. 


Bianchi et ah 


201211 . What is new here is the 


observation of the hrst UV occultation event, and the association of the BAL variability 
with reddening changes, supporting an origin of the BAL gas in the torus, a natural source 
of gas in the outer central engine of an AGN or quasar. 

Why do we observe this unusual variability in WPVS 007 and not in other BAL line 
quasars? WPVS 007 has an anom alously low luminos ity and associated black hole mass. 


size and time scales for a BALQ flLeighlv et ah 


200911 . Therefore it is possible that we 


observe behavior in WPVS 007 that is common in BALQs but cannot usnally be observed 


on human timescales. For example, in 


Leighly et al. 


(120091 1. we compare WPVS 007 with 


LBQS 1212-1-1445, a more typical BAL quasar that has a similar far-UV spectrum, but is 
100 X more luminous, with 10 x longer timescales. So, a similar occultation event in such a 
luminous object might last 600 days, while the secular decrease might last 100 years. 


While our scenario is attractive, it needs further testing. The Swift monitoring reveals 
that WPVS 007 is emerging from its occultation event. If our scenario is correct, an increase 
in flux and associated decrease in reddening would be associated with an increase in BAL 
Vmax and broadening of the C IV line as the central engine is again revealed to our line of 
sight. 
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Fig. 1.— Left: long timescale UV variability of WPVS 007 using Swift hlter UVM2. The 
upper panel shows a secular decrease since ~ 2010; the HST observations are marked by 
vertical dashed lines. The lower panel shows the variable E{B — V) derived from htting the 
photometry. The inset shows the occultation event; the vertical bars demarcate the 60-day 
minimum. Right: color-magnitude plots. The gray lines show the best ht for SMG redden¬ 


ing (solid) and intrinsic s pectral variabi 


spectrum index ay = 0.08, 


Koknbo et al. 


ity (d ashed line; best-htting power-law difference 


201J). 
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Fig. 2.— Left: the five HST spectra, plus the Swift photometry during the four recent 
observations. The spectral variability resembles variable reddening. Right: the dramatic 
variability of the C IV BAL in WPVS 007. The persistent mini-BAL near 1550 A originates 
in a separate component lying at large distances from the central engine (Leighly et ah 
2009). The maximum velocity of the broad absorption line appears inversely correlated with 
the continuum flux. The broad lines do not go below the continuum in the 1996 dim state, 
and the 2013 BALs have approximately fiat, saturated troughs, implying the presence of a 
non-variable continuum component. 
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Fig. 3.— Left panel: difference spectra illustrate the emission-line and absorption-line vari¬ 
ability; relevant lines’ restframe wavelengths are shown by vertical dashed lines. The differ¬ 
ence between the 2010 and 2013 spectra and the 2015 spectrum are normalized to zero-flux 
continuum by subtracting a constant. The broader C IV line in the earlier spectra is seen 
in the positive residuals bracketing the miniBAL (top). The Lya red wing in the low-state 
spectrum is seen in the negative residuals in the lower panel. Right panel: the spectra (iA; 
10“^^ erg s“^ cm“^A“^) (black) overlaid with the best htting model (red) and the intrinsic 
continuum plus emission lines (blue). 
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Fig. 4.— Size scales in WPVS 007 as a function of plausible black hole mass. See text for 
details. The dynamical radius of the occultor, e.g., 7?Afes(2i5o)) the dynamical radius for the 
inner BLR variability Rabs{civ), and the radius of the BAL gas based on physical properties 
Rbal are all consistent with an origin in the torus (between the inner edge, Rtk and emission 
peak Ri 2 f,m)- 
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Torus 


Fig. 5.— A simplified illustration of the scenario presented in §3.41 a) An example of a torus 
with variable scale height; realistically, the amplitude and azimuthal scale of the variation 
would be much smaller, b) An illustration of BAL winds, accelerated from the torus material 
like spray from the crest of a wave, with velocity increasing with distance beyond the torus. 
Our effective line of sight changes as material with different scale heights orbits into and 
out of our line of sight, naturally producing high velocity lines with low reddening and vice 


versa. 
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